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Abstract. New results with increased statistics arc presented for heavy flavour pro- 
duction at Q 2 ,55 150 GeV 2 and in the photoproduction limit Q 2 — ► 0. Cross sections 
for D* production, F%, the gluon density in the proton, and inelastic J/ip production 
are discussed and compared to theoretical calculations. A first measurement of the bb 
cross section is reported. 



Introduction 

At HERA positrons of 27.5 GeV collide with 820 GeV protons yielding a center 
of mass energy of 300 GeV. Heavy flavours are predominantly produced in pairs by 
photon gluon fusion. Charm quark production is expected to be a factor 200 more 
abundant than bottom quarks at this energy. Heavy flavour processes give new 
opportunities of studying perturbative QCD at center of mass energies roughly a 
factor 10 higher than in fixed target experiments. 





FIGURE 1. Charm production via photon gluon fusion (left) and J/ip production in the color 
singlet model. 



x > Contribution to "Workshop on Heavy Quarks at Fixed Target", Fermi National Accelerator 
Laboratory, Oct. 10, 1998 
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FIGURE 2. 

RAPGAP [2]. 
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ZEUS: Reconstructed D* related quantities compared to Monte Carlo simulation 
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FIGURE 3. HI: Cross sections (2 < Q 2 < 100 GeV , 0.05 < y < 0.7, p% > 1.5 GeV, 
\rj D | < 1.5) compared to NLO calculations [3] . The shaded band represents the uncertainty 



There are several methods to tag heavy flavours: "Open" charm production is 
tagged via reconstruction of D* (HI and ZEUS) or via semi-leptonic decays to 
electrons (ZEUS), b quarks have been measured via semi-muonic decays by HI. 
Finally "hidden" charm is studied via reconstruction of J/ip (see fig. 1). ip(2s) and 
T Mesons have as yet been reconstructed only in diffractive processes [la] at HERA 
and will not be reported here. 

The integrated luminosity delivered by HERA has steadily increased over the 
years. This review will cover data from 1995 (~ 6pb _1 ), 1996 (~ 10pb _1 ) and 
1997 (~ 26 pb -1 ). Most results are preliminary. 

The usual kinematic variables for deep inelastic scattering are used: 
8 = (k + P) 2 ; Q 2 = -(k - k') 2 ; x =^-- y = ^- W ' p =(q + Pf = sy - Q 2 
where k and P are the four vectors of incoming electron and proton, and q of 
the exchanged photon. 



Determination of F£ 

The inclusive cross section for production of charm in deep inelastic scattering 
(DIS) can be written as 
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FIGURE 4. F£ as function of Q 2 at fixed x (scale factors applied for clarity); F£ as function of 
x in bins of Q 2 . For comparison a NLO prediction using the GRV94HO parton densities is shown 
with the uncertainty due to the charm mass. 
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where the contribution due to Fl has been neglected since it is expected to be 
small. 

Charm is tagged through reconstruction of D* + — ► D tt + with subsequent decay 
D° — > K~n + and also the charge conjugate decay. ZEUS has presented a new 
analysis of semileptonic charm decays c — > e + + X. The electron was identified 
using the electromagnetic calorimeter and the specific energy loss dE/dx in the 
driftchamber. Details of the analysis from HI and ZEUS can be found in [lb]. 

For D* production a comparison to the RAPGAP Monte Carlo [2] simulation is 
shown in fig. 2. Reasonable agreement is found. After unfolding detector effects 
cross sections are obtained in a restricted kinematical region, examples from HI 
data are shown in fig. 3. The data are compared to a NLO calculation by Harris 
and Smith [3] using the Peterson fragmentation function. The agreement is good 
and the extrapolation to the full kinematic region is done with this calculation. 

The resulting Q 2 and x dependence of F% is shown in fig. 4. The data span Q 2 
values from 1.8 to 130 GeV 2 and 5- 10~ 5 < x < 0.02. The agreement of the different 



data sets is reasonable within errors. Also shown is the theoretical NLO calculation 
using the GRV94-HO parton density functions which reproduces the data well. A 
strong rise of _F 2 C towards low x is observed at fixed Q 2 and in Q 2 strong scaling 
violations are seen at fixed x. F£ gives a contribution of between 10% (low Q 2 ) 
and 30% (high Q 2 ) to the inclusive F 2 at an x ~ 5 ■ 10~ 4 . 

Photoproduction of D* 

When the exchanged photon is almost real, contributions due to its hadronic 
nature have to be taken into account ("resolved" processes). In NLO QCD cal- 
culations an unambiguous separation of the direct process (fig. 5a) and resolved 
processes (b and c) is no longer possible, only the sum of the two is well defined. 
There are two approaches to calculate the photoproduction cross sections in next 
to leading order: 

The "massive" approach [4] where only the light quarks u, d and s and gluons 
are active partons in the photon (and proton). Charm is only generated in the 
hard subprocess (see also fig. 5b). This approach is valid for m c ^> Aqcd- The 
"massless" approach [5,6] where also charm is an active flavour. This approach is 
valid at p t ~^> m c . 

The high statistics data from ZEUS [7] are shown in fig. 6. They are found to be 
above the massive and massless calulations. The comparison of HI data [lb] with 
massive calculations shown in fig. 7 is satisfactory. 

ZEUS has presented an analysis of D* events which contain two jets [7]. In these 
events the observed momentum fraction x° bs can be calculated, which describes the 
fraction of the photon energy contributing to the production of the two jets. A 
significant tail at low x° bs is found in the data. In the generator HERWIG this tail 
can be described by charm excitation in the photon, while considering only light 
flavours leads to discrepancies. 

Gluon density from D* events 

HI extracted the proton's gluon density function in DIS (2 < Q 2 < 
100 GeV 2 , 0.05 < y < 0.7) and in photoproduction (0.02 < y < 0.32; 0.29 < 
y <0.62; Q 2 ~ 0) [lb]. 

The observed momentum fraction x obs of the gluon is reconstructed from the 
kinematics of the final state and a differential cross section da/dx° bs is determined, 
which for the DIS data is shown in fig. 8. The correlation of x obs with the true x g 
as given by the NLO QCD calculations [3] - also shown in fig. 8 - is used in an 
iterative unfolding procedure to obtain da/dx g . The gluon density is then obtained 
by reweighting the calculation with the measured cross section. The result is shown 
in fig. 8 as a momentum distribution x g ■ g(x g ). The range 10~ 3 < x g < 0.02 is 
covered. The data from photoproduction and DIS agree well within the large errors. 
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FIGURE 5. Generic diagrams for a) direct, b) resolved charm production. In c) a NLO diagram 
is shown (flavour excitation). 
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FIGURE 6. da/dp± 
Q 2 < 1 GeV 2 and \ V D " \ 
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channels. The curves represent "massless" and "massive" calculations as indicated. 
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FIGURE 7. Photoproduction of D* from HI: da^/drj for p± > 2.5 GeV and da~> p /dp ± for 
\y D | < 1.5, where y D is the rapidity. The histograms show the NLO QCD predictions cal- 
culated according to [4] with the MRSG (solid) and MRSA' (dashed) proton parton density 
parametrizations . 
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FIGURE 8. (a) Differential cross section for ep — » eD*X as function of xf s in the visible range 
compared to the NLO QCD prediction using the CTEQ4F3 parton distribution. The shaded band 
reflects the uncertainty due to the charm mass 1.3 < m c < 1.7 GeV. (b) Correlation between x° bs 
and true momentum fraction x g . (c) Comparison of the gluon densities obtained from the two D* 
analyses in DIS and photoproduction. The systematic error is dominated by the uncertainty of 
the charm quark mass and the fragmentation parameter. For comparison the HI QCD analysis 
of the inclusive Fi measurement (shaded band) and the CTEQ4F3 parametrization are shown. 

They also agree with the result from an analysis of scaling violations in the inclusive 
measurement of Fi [8]. 

bb Production 

Due to the higher mass of the b quark the total cross section for bb production 
is expected to be 200 times smaller than that for cc production. The theoretical 
uncertainties in calculating the next to leading order predictions are, however, 
smaller [9] . HI determined the cross section for the first time in the HERA energy 
range using semi-muonic b decays [lc]. 

A photoproduction event sample was selected containing two jets of transverse 
energy > 6 GeV and a muon of transverse momentum (relative to the beam 
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FIGURE 9. Left: Definition of Vt re i- Right: The measured rel distributions and the result 
of the fit (solid line) ; the contributions from beauty, charm and background are shown separately. 



direction) p£ > 2 GeV in the central detector region 35° < 9^ < 130°. 

The thrust axis 2 was determined for each jet in order to approximate the b flight 
direction. The transverse momentum of the muon p^ reX with respect to the jet 
is used as a discriminating variable: Muons from b decay show a Pt rei spectrum 
extending to higher values than c decays (see fig. 9 for an illustration of the method). 

The background comes from the production of the light quarks u, d and s, which 
is roughly a factor 2000 larger than b production. Punch through and decay in 
flight lead to false muon signatures. The contribution is determined from data 
using an independent dataset and using the muon fake probability and the hadron 
composition from a well tuned and checked simulation program. The resulting p^ rd 
spectrum is shown in fig. 9 indicating the background contribution (23.5%), which 
is absolutely determined. The fractions of b and c quarks are obtained from a fit 
to the data distribution, yielding (51.4 ± 4.4)% and (23.5 ± 4.3)%, respectively. 

The cross section in the visible kinematic ranee of Q 2 < 1 GeV 2 ; p^ > 2 GeV; 
95 < W 1P < 270 GeV; 35° < 0** < 130° is determined as 

a(ep ^ebb + X) vis = 0.93 ± 0.08±g; 2 £ nb, 

where the first error is statistical and the second systematic. Contributions to the 
systematic errors are the branching ratio b — * Xfiv, the energy scale of calorimeters 
and detector efficiencies. 

2 ) The thrust axis is the axis which maximizes T — max(|^), where the sum runs over all 

particles belonging to the jet and pf is the component of the particle momentum parallel to the 
thrust axis. 
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FIGURE 10. Inelastic J/t/j photoproduction: a) cross section as function of W lp for 
0.4 < z < 0.9 and px > 1 GeV. The lines correspond to the NLO predictions from [13] with 
GRV or MRSA' [20] parton density functions (m c = 1.4 GeV and A QC d = 300 MeV). b) da/dz 
for 50 < W 1P < 180 GeV and pr > 1 GeV. The NLO computation is shown as a solid line. The 
dashed and dotted lines are given by the sum of the colour-singlet and the colour-octet leading 
order calculations performed in [15] and in [16]. The latter include estimates of higher order QCD 
corrections due to initial state radiation. 



The corresponding direct LO cross section from the AROMA [10] simulation is 
0.19 nb, roughly a factor 5 lower. The fraction of c quarks determined from this 
analysis leads to the same cross section for ep — > eccX as previously determined 
from analysis of D* production [11]. 



Inelastic J/ift production 

New data on charmonium (J/ip, ip') and T production have been presented by 
HI and ZEUS [la]. Here we will concentrate on "inelastic" J/ip production as 
opposed to the diffractive processes which dominate the cross section at low Q 2 . 
Inelastic J/ip production could at lower W 1V (fixed target regime) be well described 
by the Colour Singlet Model (CSM). For HERA the CSM cross section calculations 
are available in NLO [13]. 

As is well known the CSM fails to describe charmonium production in pp collisions 
at high pt [12]. Colour octet contributions have been proposed for an adequate 
description. The NRQCD factorisation approach (NRQCD = Non Relativistic 
QCD) describes any process A + B^J/tp + X&s& sum over colour singlet and 
colour octet contributions. 

Whereas the transformation of a colour singlet 3 Si state into a J/ip can be calcu- 
lated using the measured semileptonic decay width, the transition of a colour octet 
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FIGURE 11. Differential cross sections for the inclusive (open points) and inelastic 
(Ef w d > 5 GeV, black points) ep — ► e J/i/> process, a) da/dQ 2 , b) da/dp 2 c) da/dz and 
d) da/dy* . The curves are predictions [17] within the NRQCD factorization approach for the 
colour singlet contribution (dashed) and the sum of singlet and octet contributions (full line). 

state to J /if) is non-perturbative and at present not calculable. Therefore predic- 
tions for the cross section at HERA use the non perturbative transition matrix 
elements extracted from the CDF data. 

The ZEUS collaboration has updated their photoproduction data [Id]. The re- 
sults for the 7p cross section as function of W lp and of z are shown in fig. 10. The 
data agree well with the next to leading order pQCD calculation in the color singlet 
model [13]. The variable z is defined as z = ^jf^- m where the latter approxi- 
mation holds in the proton rest frame. In fig. 10 b in addition to the CSM in NLO 
calulations using the NRQCD/factorisation approach [14-16] are shown. The up- 



per curve was calculated in LO using the transition matrix elements extracted from 
CDF data in LO and shows a strong rise towards high z values. The lower curves 
also take into account higher orders approximately as explained in refs. [14-16]. 
Doing so leads to modifications in the non perturbative matrix elements and/or in 
the cross sections themselves. The net effect is a decrease of the predicted rise at 
high z. 

HI has for the first time determined the cross sections for inelastic J/ip pro- 
duction at Q 2 > 2 GeV 2 [Id]. The results are shown in fig. 11. Two data sets 
are shown, a completely inclusive one (open points) and one where the diffractive 
contributions have been removed by a cut on the energy in the forward region of 
the detector as suggested by Fleming and Mehen [17], whose LO calculations are 
shown for comparison. The data are seen to be far above the CS contributions. 
The magnitude of the data is reproduced better taking into account colour octet 
contributions. The shape of the latter leaves, however, much room for improve- 
ment, in particular in the rapidity y* in the 7*p center of mass system. Note that 
the NRQCD calculations are performed at the parton level, no smearing due to the 
transition into J/ip is taken into account. 

Summary 

Due to increased statistics detailed analyses of heavy flavour production in ep 
collisions are performed in a variety of channels and kinematic regions, bb pro- 
duction was observed for the first time in photoproduction via semi-muonic decay 
of the b-quark. The cross section was found to be considerably larger than the 
leading order predictions for the direct process. 

In the range 2 < Q 2 < 130 GeV 2 cross sections and the charm contribution to F| are 
determined and found to agree with next to leading order predictions. In photo- 
production the validity of different approaches to calculate next to leading order 
corrections is being studied in various kinematic regions. 

Since photon gluon fusion is the dominant process a direct determination of the 
gluon density in the proton was carried out in DIS and in the photoproduction 
regime. The result agrees with the indirect determinations from scaling violations. 

Inelastic J/ip production is studied in photoproduction and DIS and is well de- 
scribed in photoproduction by the color singlet model alone in next to leading order. 
In DIS the data have been compared to LO color singlet and color octet predictions 
(at parton level). In the latter rough agreement in absolute normalisation is found, 
while the color singlet model reproduces the shape of the data slightly better. 

Acknowledgement I wish to thank the organisers for a very pleasant and fruit- 
ful meeting and my colleagues at ZEUS and HI for supplying their data and for 
discussions. 
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